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The status and ongoing developments in the field of deep inelastic scattering pre- 
sented at the DISOl workshop in Bologna are discussed from both the experimental 
and the theoretical perspective. 

1 Introduction 

In the last year considerable progress has been made in many aspects of the 
study of the structure functions of photons and nuclei. The most important 
results presented in the structure function session at the DISOl workshop in 
Bologna are reviewed. The new theoretical developments are discussed in 
Section ^j. New experimental information on photon structure is summarised 
in Section ^ while Section ^ reviews the new experimental results pertaining 
to the structure of the proton. The material presented here is intended only 
as an overview, and reflects the personal view of the working group convenors. 
Many of the important details of the investigations have of necessity had to be 
omitted. The reader is referred to the write-ups of the individual presentations 
that are to be found elsewhere in these proceedings. 

2 Theory 

It is remarkable that a theoretical framework that was established more than 
twenty years ago - short-distance factorisation and next-to-leading (NLO) 
DGLAP - still gives a very good description of almost all deep inelastic struc- 
ture function data, for example see Section ^ below. In recent years attention 
has focused on calculations that attempt to improve the theoretical predic- 
tion by going beyond the standard 'massless quark NLO-DGLAP' framework. 



*INVITED TALK PRESENTED AT THE 9TH INTERNATIONAL WORKSHOP ON 
DEEP INELASTIC SCATTERING (DIS 2001), BOLOGNA, ITALY, 27 APR - 1 MAY 
2001, TO APPEAR IN THE PROCEEDINGS. 
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Previous DIS Workshops reported the significant progress on incorporating 
non-zero quark masses (particularly rric) into the analysis, and the resulting 
improved calculations are now used routinely in the experimental analyses. 
The outstanding theoretical deficiencies and uncertainties can be classified 
into various types. Within leading- twist perturbation theory, the obvious goal 
is a complete NNLO description. At least for massless quarks, this seems to be 
within reach (see below) and the likely phenomenological impact has already 
been investigated. A variation of this approach is to focus on the kinematic 
a; — > 0, 1 limits and resum the leading logarithmic corrections to all orders 
in perturbation theory, to investigate whether the fit to data is improved. 
Non-perturbative corrections arise in various ways, for example higher twist 
(l/Q^) , n > 1 power corrections, and also heavy nuclear target corrections. 
A number of theoretical contributions to this session addressed these issues, 
and the main results of these studies are summarised briefly below. 

Parton distribution functions (PDFs) extracted from deep inelastic scat- 
tering structure functions data play a central role in the calculation of hard 
scattering cross sections at the Tevatron and LHC. A precise knowledge of 
the PDFs is absolutely vital for reliable predictions for signal and background 
cross sections. In many cases, it is the uncertainty in the input PDFs that 
dominates the theoretical error on the prediction. Such uncertainties can arise 
both from the starting distributions themselves, reflecting the uncertainties 
in the data input to the global fit, and from evolution to the higher scales 
typical of hadron collider hard scattering processes, which is sensitive to un- 
certainties in as, unknown higher-order corrections, other parton flavours, etc. 
By way of illustration. Figure || shows the (a;, Q^) values corresponding to the 
production of heavy objects (e.g. a VF or Higgs boson, a ti pair, a multijet 
final state etc.) of mass M and rapidity y. We assume leading order kine- 
matics, so that X — M exp{zLy) / and Q ~ M. As an example, a W boson 
(M = 80 GeV) produced at rapidity y — 3 corresponds to the annihilation of 
quarks with x = 0.00028 and 0.11, probed at = 6400 GeV^. Notice that 
in this example, quarks with these x values are already 'measured' in deep 
inelastic scattering (at HERA and in fixed-target experiments respectively), 
but at much lower Q^. 

A rigorous and global treatment of PDF uncertainties remains elusive, 
but there has been a significant advance in the last few years, with several 
groups now introducing sophisticated statistical analyses into global or quasi- 
global fits. Further progress was reported at this meeting and is summarised 
below. 
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2.1 Higher orders, higher twists 

In order to match the precision of present and future DIS data, it is necessary 
to go beyond NLO in DGLAP perturbation theory. The most relevant quan- 
tities here are the NNLO contribution to the DGLAP sphtting functions, i.e. 
the functions P'^^ in the expansion 

Sphtting Function: P(.t, a^) = P^"^ + a^P^^' (x) + a^P^^^ (x) + ... (1) 

together with the corresponding coefficient functions: 



Coefficient Function: 



(2) 
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The NNLO contributions to the latter are easier to calculate and are already 
known for many of the cross sections and structure functions of interest.^ 

Using recent iV = ID, 12 moment calculations by Retey and Vermasereni, 
van Neerven and Vogtcl have updated their approximations for P^'^^{x), in 
which the moment information is combined with sum rule and symmetry con- 
straints and known leading behaviour for a; ^ 0, 1. The resulting approximate 
functions are certainly sufficient for x > 10~^ phenomenology, and probably 
adequate for x > 10~^ as well. Further_information can be found in the 
contribution by Vogt to these proceedinggj. It was reported by Vogt at the 
meeting that the exact calculation of the NNLO splitting functions should be 
completed by the end of the year. 

Fully quantitative predictions for higher-twist contributions to DIS struc- 
ture functions remain elusive. A potentially important twist-4 contribution 
comes from the 4-gluon operator matrix ^filement combined with the 4 gluon 
— > 2 quark coefficient function. Bliimleino presented an estimate of this con- 
tribution to the structure function slope 9^2/9 log Q^. The calculationd is 
performed using time-ordered perturbation theory and is limited to dp^/p^ 
accuracy in the calculation of the corresponding Feynman diagrams. Numer- 
ical results are obtained by modeling the 4(7 distribution amplitude. At small 
X it is the screening (rather than anti-screening) terms which dominate and 
reduce the increase of the slope due to the leading twist-2 contributions. The 
numerical results, see Figure |[ show that this effect is significant only below 
cc ~ 2 • 10~^ for Q"^ ^ b GeV^. This is beyond the kinematic range that can 
be probed at HERA but may be accessible at future lepton-hadron colliders 
operating at higher energies. 

Another type of power correction, of more universal origin, can arise from 
analytic resummation of perturbative contributions to all orders. Such resum- 
mations typically lead to expressions of the form 

^{k^Tasie) (3) 

which are ill-defined because of the Landau pole in the running coupling. Ex- 
panding cks(fc^) in powers of as(Cp) yields a perturbation series for / whose 
coefficients grow as n\. The resulting ambiguity in the resummed perturbative 
result signals a nonperturbative power correction of the form (A^/Q^)". An 
alternative way to tackle this problem, in the context of QCD quark form 
factors, using dimensional regularisation was presented by Magneafl (see also 

"Note however that a^"^^ is still not calculated for some of the key cross sections used in 
global fits, e.g. da/dEi^^ , da/dy^y, d'^a^^ /dM dy, etc. 
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0) . If the calculation is performed consistently in d = 4 — 2e > 4 dimensions 
then the pole moves off the real axis and analytic resummed expressions for 
the form factors can be obtained in terms of gauge and renormalisation-group 
invariant quantities. Applying these techniques to the physically relevant DIS 
structure function (N) moments yields the standard expression cx N{h? /Q'^) 
for the leading power correction at large N . The generalisation to more com- 
plicated QCD amplitudes, for example those involving at least two coloured 
particles, is more difficult and more work needs to be done. 

Small- a; structure function phenomenology incorporatiag higher-twist 
contributions was the subject of a presentation by KotikovEl. The starting 
point of the analysis is the set of analytic solutions of^lhe DGLAP equations 
at small x obtained from a flat parton input at Qq lI. These solutions are 
supplemented by higher-twist (twists 4 and 6) contributions as obtained in a 
renormalon modetj. A very reasonable few-parameter fit to HERA structure 
function data is obtained (see 13), suggesting that higher-twist contributions 
do contribute to the observed rise of the F2 structure function at low values 
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of X and Q^. l— . 

Another example of QCD resummation was discussed by ErmolaevO 
in the context of the small- a; asymptotics of the non-singlet unpolarised 
and polarised structure functions. Accounting for both double- and single- 
logarithmic contributions and using a fixed QCD coupling yields the asymp- 
totic behaviour, 



^r.ffP-ir^/^i (4) 




a+ = V2asCf/7r, a_ = a+ ^1 + 1/{2N^) (5) 



The situation is more complicated when oip-is allowed to run with the internal 
ladder transverse momenta as argumentalS. Two-dimensional infrared evo- 
lution equations are constructed and solved numerically. The exponents a± 
become functions of Aqcd and fj,, the infrared cut-off parameter. Quantita- 
tively, a± ^±{i^) = 0.37/0.4 at 1 GeV, which in fact is in good agreement 
with the phenomenological values extracted from fits to xF^ data at small x. 

Finally, it is by now well established that it is difficult to disentangle 
the underlying small-x QCD dynamics (BFKL, CCFM, DGLAP, etc.) from 
inclusive quantities (e.g. F2) alone. More exclusive quantities, for example 
transverse energy flow and dijet productionxates and correlations, offer a mate 
decisive test, at least in principle. SzczurelJl3 presented the results of a studyllj 
of how dijet azimuthal correlations da/d(j)jj in real and virtual photoproduc- 
tion at HERA can be used to probe the unintegrated gluon g{x,k'j_) in the 
proton. The benchmark distribution used iiLthe study is the two-component 
(soft -I- hard) model of Ivanov and Nikolaevtj, which leads to distinctive kine- 
matical dependence of the dijet correlations. It will be important to make 
comparisons with the predictions of other models, and of course with data. 

2.2 PDF fits and uncertainties 

Progress in obtaining an improved quantitative understanding of parton dis- 
tribution functions and their uncertainlies was reported at the meeting. In 
the context of 'best global fits', ThorneE^ reported on a recent update of the 
MRST (NLO-DGLAP) global fit. The main new ingredients in the fit are (a) 
new small- a; structure function data from HERA, and (b) new high-i?T jet 



data from the Tcvatron (see Section 4.2 below for more details). For the latter, 



a proper treatment of the correlations between the various systematic errors 
is essential and is incorporated in the MRST fitting procedure. The main 
effect of the new data is to further constrain the gluon distribution at both 
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small and large x. As for previous global analyses, the overall fit is very good, 
but there are some interesting new features: (i) the negativity of the small- a; 
(NLO, MS) gluon near the starting scale at 1 GeV^ is now firmly established 
(although the gluon becomes positive for all x > 10^^ for > 2 — 3 GeV^); 
(ii) there is a slight 'tension' between the as(Mf ) values preferred by the DIS 
(0.121) and jet (0.117) data. p 

Turnuifi_to uncertainties on (or due to) PDFs, the MRST approachO 
(see also tM3) has been to focus on the uncertainties on particular physical 
quantities Saohs (for example, the Higgs or weak boson production jCrpas-sec- 
tions), rather than on uncertainties on the PDFs themselves, Sfi l23'EllE2lBE3. 
Thorne reported on a new determination of the PDF uncertainty on aw at 
the Tevatron and LHC. A priEji'ious investigation of this uncertainty using the 
Lagrange Multipher methodES (see below) had resulted in ^ctpdf — 5 — 8%. 
The MRST approach is to determine which parts of the global fit fail when 
trying to force higher/lower (Jobs, and to impose the criterion that no individ- 
ual data set has a less than 1% confidence level. This leads to an estimate 
of the PDF dependent uncertainty in the total W cross section of ~ ±2% 
at both the Tevatron and LHC. An example of how one of the components 
of the MRST global fit fails when aw is forced up and down from its central 
prediction is shown in Figure ||. 

Recent work on parton distributionjmcertainties in the context of the 
CTEQ global fit was reviewed by StumpEa. Two methods iiave been adopted. 
In the Lagrange Multiplier constrained fitting approachcS the global fit is 
repeated for successive constrained values of a physical observable trx, and 
a Xgioijai "vs- ax profile is established. The issue is then what 'tolerance' 
(increase in Xgiobai from the minimum) corresponds to a standard confidence 
level. In the CTEQ study, the tolerance is determined by considering the 
response of the individual experiments to the sample PDFs. The 90% CL 
allowed ranges for each experiment are then combined into a single overall 
tolerance. The second 'Hessian Matrix' approach follows a more traditional 
error analysis by focusing on the allowed ranges of the individual parameters 
ai, ...ad that define the PDFs: 

d 

Xgiobai ~ - a°) + • ■ ■ (6) 

Assuming a quadratic approximation, the error on a given observable X is 
then 



7 



CCFRXF, data . W^^(\.02\). W^y(0.973), W-jgy(1.034), W^y(0.956) fit 




1.2 - 




K-0.75 (+0.1) 



1 10 10" 

Q- (GeV^) 

Figure 3. Comparison of CCFR i^3(x,Q'^) data with theory for (Tv(/(TEV) changing in 
magnitudej-fejy factors of 1.021, 0.973, 1.034 and 0.956, from the new MRST study reported 
by Thornetfl. 



The eigenvectors of the Hessian matrix H define orthogonal directions in 
parameter space that can be used to define a set of {2d + 1) PDF basis sets. 
Again there is the issue of the definition of the standard tolerance, i.e. the 
allowed variation of Ax^ in Equation |^. The PDF basis sets can also be used 
to define bands of uncertainty for individual PDFs as the envelope of the 
extreme curves at given values of x. An example is shown in Figure ^. Notice 
how the band shrinks as increases, and also how the choice of parametric 
form causes a slight (artificial) decrease in the band width between the low 
and high x regions. 

Finally, in a thought-provoking contribution Collin£3 (see also EZl) ad- 
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Figure 4. Extreme glpyn distributions at Q = 2 and 10 GeV, from the CTEQ Hessian 
Matrix error analysiscj- The curves show the extremes at some particular x, while the 
shaded regions represent the envelope of such extremes. 



dressed the issue of how to determine how 'good' a good global fit really is. 
The point is that a reasonable overall Xtot '^^^ mask the fact that certain 
subsets of data (i.e. individual experiments) are in strong disagreement with 
the theory (which could of course have various explanations, ranging from 
errors or deficiencies in the experimental analysis or theoretical calculation to 
an unexpected contribution from new physics). The procedure is therefore to 
consider the quality of fit to individual data subsets, explore the parameter 
space to find PDF sets that are 'best fits' to these individual data sets and 
to examine by how much the fit to a particular data set is improved relative 
to the increase in overall Xtot- ^ substantial decrease in an individual is a 
symptom of a bad fit to that particular data set, which could then be further 
studied to investigate the source of the disagreement, ft is interesting to note 
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that both the CTEQ and MRST global fits contain 'problem' data subsets 
according to this definition. 



2.3 Parton distributions in heavy nuclei 




When deep inelastic neutrino scattering structure function data are used 
in a global PDF fit, account must be taken of the fact that the PDFs in heavy 
nuclei are different from those in a nucleon. It is important to quantify these 
differences, not only for use in global PDF fits but also, for example, when 
interpreting possible quark-gluon plasma signals from heavy ion collisions. 
There haizE_bem several recent analyses to determine the pajcton distributions 
in nuclein§n3'L3. KumanoES reported the results of a studyE2l based on fits to 
electron and muon F2 data on nuclear targets. A simple ansatz is used to 
parameterise the nuclear modification of the PDFs: 

w,{x,A) = 1 + {1-A-^Y"'k{x) (8) 

where hi is taken to be either a quadratic or cubic polynomial in x divided by 
a power of {1 — x). The parameters in hi are determined by a minimisation 
procedure. A reasonable overall fit is obtained (x^/dof = 580/302), with the 
excess mainly due to apparent inconsistencies between data sets, e.g. E665 
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vs. NMC. The resulting weight functions for calcium are shown in Figure |[ As 
only structure function data are used in the fit, there are large uncertainties 
on the extracted gluon and sea-quark weight functions. Comparisons with 
Drell-Yan, prompt-7 and heavy qnfn'V production data are needed to further 
constrain g"^, q-^, see for exampla^Jfi 



3 Photon Structure 



The investigation of the structure of the photon is a very active field of re- 
search at LEP as well as at HERA, seecJ for a recent review. At this confer- 
ence new results from LEP and HERA have been presented which are briefly 
summarised here. 

New investigations on the structure £if the photon based on jet pmduc- 
tion were presented by HI for quasi-realo, and by ZEUS for virtuaO pho- 
tons. The new HI result, Figure |[ is consistent with the predictions based 
on existing parametrisations of the hadronic photon structure function F2 
and at present the data are not precise enough to distinguish between differ- 
ent parametrisations. This has to be confronted with the earlier result from 
ZEUSE3 which suggested that the parametrisations of , obtained from fits 
to measurements made at e"'"e~ colliders, are too low for medium values of 
Bjorken x and at factorisation scales of several hundred GeV^. There are sev- 
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Figure 6. Measurement of the structure of quasi-real photons from HI. 



11 



eral differences between the ZEUS and HI analysis such as the choice made 
for the value of as, the parton distributions used for the photon, and most 
notably the difference in the corrections applied to the data. The HI data 
are corrected for detector as well as hadronisation effects and are shown at 
the partonic level. In contrast, the ZEUS results are corrected only for de- 
tector effects and phase space regions are selected, where the hadronisation 
corrections, as implemented in Monte Carlo models, are found to be small. It 
remains to be seen how much of the apparent differences between the results 
can be explained by the different analysis methods. 

Also for the measurement of the structure of virtual photons some clari- 
fication is needed. The recent ZEUS measurement, Figure 0, indicates some 
shortcomings of the prediction from the SaScZI parametrisations. The sup- 

ZEUS PRELIMINARY 



•ZEUS 96/97 yp 
— SaS1D(HERWIG 5.9) 
■■■ GRVLO(HERWIG 5.9) 




Figure 7. Measurement of the structure of virtual photons from ZEUS. 

pression of the photon structure with the photon virtuality is studied based 
on the ratio of the cross sections for low and large values of x. The SaSlD 
prediction fails to describe this ratio. Again, there seems to be a difference 
between the ZEUS result and the earlier result from HI on the structure of 
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virtual photonsE^I. The HI result, which is valid for an almost identical phase 
space region, showed good agreement between the data and the predictions 
based on the GRV parametrisation of F2 with a virtuaUty suppression as 
given by the Drees-Godbole scheme. 

It is certainly desirable to complement the measurements of F2 with the 
jet measurements from HERA which extend to larger factorisation scales, 
when fits for the parton distribution functions of the photon are performed. 
However, first it has to be seen if a consistent picture of the various HERA 
results can be established. 
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Figure 8. Prospects for the measurement of the at a future Unear coUider, from 
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Recently at LEP progress has been made in the measurements of 
{x,fB). The phase space of the measurement has been extended by 
both to lower values of x and to larger values of . In addi- 



tion, significant progress in reducing the systematic uncertainty has been 
achieved by using improved Monte Carlo models to describe the hadronic 
final statecj, as well as by utilising a more sophisticated unfolding procedure 
in two dimensionscll. In view of these improvements a re-analysis of all LEP 



13 



data using the newly available analysis tools would significantly improve on 
the precision of F2 ■ 

The prospects of future investigations of the photon structuxe in the con- 
text of the planned linear collider programme are very promisingfj. The e+e~ 
linear collider will extend the available phase space, as shown e.g. in Figure ^ 
for the measurement of the evolution of F2 at medium x. The higher 
beam energy and luminosity compared to LEP also allows for the investiga- 
tion of novel features like the measurement oLthe flavour dependence of -F^ 
by exploring the exchange of Z and W bosonsa. When using highly polarised 
beams tlie measurement of structure functions can be extended to polarised 
photonaU. The investigation of the photon structure will enter a completely 
new level of sophistication if a photon collider can be realiseclEJ. Then photons 
of large energy with a rather moderate energy spread could be brought into 
collision, instead of the presently used soft photons from the Bremsstrahlungs 
spectrum at electron-positron colliders. 



4 Proton structure 

4.1 Deep inelastic scattering at high-Q^ 

New measurements of the high-Q^ neutral current (NC) and charged current 
(CC) deep_^inelastiC|_scattering (DISl^^'^p cross sections were presented by 
the ZEUSB and HIEI collaborations 113. Figure | shows the differential cross 
section da/dQ^ for NC DIS. The ZEUS e+p data, presented for the first time 
at this conference, were obtained using approximately 60 pb~^ of e^p data 
collected at a centre-of-mass energy of y/s — 318 GeV over the years 1999 
and 2000. HI has combined 46 pb^^ of e+p data (from a total of 65 pb~^ 
collected in 1999 and 2000) at = 318 GeV with the 38 pb"^ of data 
collected over the years 1994 to 1997 at ^ = 300 GeV. The NC cross section 
is observed to fall by roughly 6 orders of magnitude as goes from ~ 
200 GeV^ to g2 - 10 000 GeV^. For < 1000 GeV^ the NC e+p cross 
section is observed to be approximately equal to the e^p cross section. In 
this kinematic domain NC e^p DIS is mediated by single photon exchange 
and the dependence of da/dQ^ reflects the l/Q^ dependence of the photon 
propagator. As approaches the Z-boson mass squared {Q^ ~ -^I) the e^p 
cross section falls below that for e~p NC DIS. At such large both photon- 
and Z-exchange contributions must be taken into account. The interference 
between the two contributions suppresses the e^p cross section and enhances 
the e~p cross section. The StandariLModel (SM), which incorporates all these 
effects, evaluated with the CTEQ5E3 parton density functions (PDFs) gives a 
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good description of the data. 




Figure 9. The high-Q^ e^p neutral current and charged current deep inelastic scattering 
cross sections measured by the ZEUS and the HI collaborations (points). The predictions 
of the Standard Model evaluated with the CTEQ5D parton density functions are shown as 
the lines. 

The differential cross section da/dQ^ for CC e^p DIS is shown in Figure 
^. The cross sections fah slowly for 200 < < 1 000 GeV^ and then more 
rapidly as approaches the M^-boson mass squared (Q^ ~ ^w)- 
SM, evaluated with the CTEQ5 PDFs, gives a good description of the data. 
The e~p CC DIS cross section is always larger than the e+p one. The CC 
interaction distinguishes between quarks of different flavour. Thus, the e^p 
CC process picks out the up-type quarks and the down-type anti-quarks while 
e^p CC DIS picks out the down-type quarks and the up-type anti-quarks. 
The e~p cross section is larger than the e"*"p cross section because there are 
more u-valence quarks than c?-valence quarks in the proton and because the 
chiral nature of the CC interaction suppresses the contribution of the down- 
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type quarks in e+p CC DIS. Note that for > 10 000 GeV^ the four cross 
sections for NC and CC DIS in e+p and e~p colhsions all have a similar 
magnitude. This is a striking experimental verification of the electroweak 
unification embodied in the Standard Model. 

The reduced cross section (t'^q is related to the double differential e^p 
NC DIS cross section by 



'NC 



dxdQ'^ 



(9) 



where Y± — 1 ± (1 — y)^ and y is related to and x by = xys. Figure 
shows CT^J' plotted as a function of for several values of x. For < 
1000 GeV^ and a; ~ 0.1 the reduced cross section is observed to be almost 
independent of while for < 1 000 GeV^ and x > 0.25 it falls slowly 
with Q^. For above 1 000 GeV^ the e~p reduced cross section exceeds the 
e+p reduced cross section. The Standard Model (SM) expectations for (J^q' 
obtained by evaluating Equation ^ using standard parton density functions 
gives a good description of the data. The SM accounts for the difference 
between a'^Q and in terms of the structure function xF^'^ so that xF^'^ 



Both ZEUS and 



may be extracted from the difference between ct^c' ^^^d o'nc'- 
HI have extracted xF^'~^ in this way, the results are plotted in Figure [iO^D. 
A good description of the data is obtained when xF^'~^ is evaluated with 
standard PDFs. More e^p data is required to improve the precision of xF^'-^. 
With sufficient data xF^'~^ will be used to constrain the valence quark PDFs. 



The charged current can also be used to study the partonic content of the 



proton. The reduced cross section for e+p CC DIS, , is defined to be 



'CC 



2'KX 



dxdQ^ 



T Y^xF^ 



(10) 

Figure 11 shows function of x for values of between 200 GeV^ and 

20 000 GeV^ . The SM expectation of Equation [l^ gives a good description of 
the data. At leading order in QCD is given by 



CTcc^ = x[u + c+{l-yY{d + s)] 



(11) 



while is given by 



'CC 



CTcc^ ^x[u + c+{l- yf{d + s)] 



(12) 
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Figure 10. (a) The reduced cross section for e^p neutral current deep inelastic scattering 
measured by the ZEUS and the HI collaborations (points). The predictions of the Standard 
Model evaluated with the CTEQ5D parton density functions are shown as the lines, (b) 
The neutral current structure function xF^'-^ extracted from e^p neutral current deep 
inelastic scattering data by the ZEUS and HI collaborations (points). The predictions of 
the Standard Model evaluated with the CTEQ5D parton density functions are shown as 
the solid lines. 



Figure [ill shows how the different quark flavours contribute to the reduced 
cross section. It can be seen that, at high-x, (t^q is sensitive to the c?-quark 

density while o'^q is sensitive to the u-quark density. 

These results, important in themselves, represent a significant milestone 
in the development of the HERA physics programme. DISOl took place at 
a time when major upgrades to the HERA machine and experiments were 
nearing completioiO^EiJ. The HERA machine upgrade will provide polarised 
electron and positron beams for the colliding beam experiments ZEUS and 
HI at a luminosity five times greater than has been achieved to date. The 
new high-Q^ cross section results, therefore, represent the culmination of the 
programme of measurements of the unpolarised DIS cross sections at high-Q^. 

4-2 Determination of as and extraction of the gluon density 

A highlight of the structure function session was the variety of new, high 
precision, simultaneous determinations jXif the strxuig coupling constant, as, 
and the gluon density, xG. The ZEUSEJ and H1E3 collaborations presented 
new analyses in which the scaling violations of the structure function -F^^ 
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Figure 11. The reduced cross section for e+p (a) and e~p (b) charged current deep inelastic 
scattering measured by the ZEUS and the HI collaborations (points). The predictions of 
the Standard Model evaluated with the CTEQ5D parton density functions are shown as 
the solid lines. Quark contributions, evaluated using the CTEQ5D PDFs, are shown as the 
dashed lines. 



are exploited to determine xG and as ■ Both collaborations have performed a 
NLO QCD fit using their own data on o'^q together with structure function 
data from other experiments. The analyses proceed hji, parameterising the 
PDFs at a starting scale Qq. The DGLAP equationsE^ are used to evolve 
from Qq to the value of that corresponds to a particular data point. The 
relevant measured quantity is then determined from the evolved PDFs and 
the contribution to the total calculated. The ZEUS and HI analyses differ 
both in philosophy and in the details of the fitting technique used. 

The goal of the HI collaboration was to determine xG and as as precisely 
as possible uaiag as few data sets as possible. The fit includes data from Hla 
and BCDMSE3 (proton target data) in order that the lever arm in x and 
be sufficient to constrain the parameterisation. The PDFs are parameterised 
at Qq = A GeV^. The PDF parameterisation is based on effective valence 
and sea distributions together with a parameterisation for xG. The number 
of parameters in the gluon and sea distributions are chosen to saturate the 
X^- The c- and 6-quark contributions to the cross sections are evahiated using 
photon-gluon fusion using a massive fixed flavour number schema^. 

ZEUS parameterised the u- and c?-valence, the sea and the gluon PDFs at 
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Qq = 7 GeV^. In addition, the difference between the.u- and d-quark density 
{x{u—d)) was parajjaeterisedi-J^he recpit ZEUS daiaEll was used together with 
data from BCDMSy , CCFREa, NMCEI and EGfiSB A modified massive vari- 
able flavour number scheme (RT-VFN schemeEil) was used. This treatment 
of the charm quark smoothly interpolates between the low-Q^ region where 
threshold effects are important to the high-Q^ region where the charm quark 
mass is negligible. 

Both ZEUS and HI have incorporated a very careful treatment of the 
experimental systematic errors which includes the point-to-point correlations. 
HI allowed the fit to determine the best values for the various systematic 
uncertainties, the ZEUS fit did not. The resulting gluon distributions are 
compared in Figure [l^. When the differences between the two analyses are 
taken into account the agreement between the two results is reasonable. The 
corresponding values of as are 

as = 0.1172 ± 0.0008(stat. + uncor. sys.) ± 0.0054(cor.sys.) (ZEUS) (13) 

as = 0.1150 ± 0.0017(expt.)+°:°°°5 (model) ± 0.005(scale) (HI) (14) 

ZEUS breaks up the total experimental uncertainty into the contribu- 
tion arising from the correlated experimental systematic uncertainties (cor. 
sys.) and the combined statistical and uncorrelated systematic uncertain- 
ties (stat. -|- uncor. sys.). HI has combined the statistical and experimental 
systematic contributions into a single error labeled (expt.). The HI 'model' 
error accounts for uncertainties in the theoretical model such as those arising 
from the choice of the functional form or the treatment of heavy quarks while 
the 'scale' error is evaluated by varying the factorisation and renormalisation 
scales. The values obtained are consistenL_with one another and with the 
values obtained in various global analysesE3. Given the importance of the 
measurements it will be important to understand in detail the range of appli- 
cability of the. results pid, if possible, to converge upon a common approach. 

The ZEUSE^I and H1E3 collaborations have exploited the sensitivity of DIS 
jet production cross sections to determine as and xG. With ag set equal to 
the world average, xG may be determined in a fit to the measured inclusive 
jet cross_sections. This has been done by HI and the result is shown in Fig- 
ure P^E£|. Alternatively, the PDFs can be taken from one of the standard 
parameterisations and the jet cross sections used to determine ag. When ag 
is determined in this way it is essential to estiniate the joiicertainty in ag 
which arises from the choice of PDF. Both ZEUSE3 and HlEj have performed 
such analyses. The HI analysis uses the inclusive jet cross section. ZEUS 
has chosen to analyse the ratio of the cross section for the production of 2- 
jets (plus the target remnant) to the total inclusive DIS cross section. The 
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Hl+ZEDS 




Figure 12. (a) The gluon distributions extracted from the ZEUS and HI NLO QCD fits 
to structure function data. The statistical and systematic error bands are shown shaded. 
Details of the extraction of these error bands may be found elsewhere in these proceedings, 
(b) The gluon distribution extracted from a NLO QCD fit to the HI inclusive jet cross 
section data. 



use of the ratio reduces the sensitivity of the result to the choice of PDF. A 
compilation of the values of ag obtained from jet production data is shown in 
Figure |l^ together with the results obtained from the scaling violations of F2. 
The results agree with one another and with the world average. The way in 
which as depends on Q^j-the running of as, may also be investigated using 
the jet data. Both ZEUa^ and H1E3 have performed such analyses and the 
results are consistent with the QCD expectation. Finally, it is also possible 
to estimate xG and as simultaneously from the jet data in a fit conceptually 
similar to that jjsed in the analysis of the scaling violations of F^'^ described 
above. The H1E3 collaboration has performed such an analysis. The fit results 
in estimates of the gluon distribution and as which are correlated. 

The techniques discussed above do not constrain the gluon distribution 
for a; > 0.1. A promising new development is the measurement of the double 
differential inclusive jet cross section (Pa / dExdrj (where Et is the jet txana- 
verse energy and 77 is the jet pseudo-rapidity) by the DO coUaboratiorEJ'Ea. 
The cross sections, shown in Figure have been measured over the kinematic 
range |?7| < 3 and 50 < Et < 500 GeV. Since the beam energy at the Tevatron 
is 900 GeV, these cross sections are sensitive to xG for values of x as high as 
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Figure 13. Compilation of the various determinations of the strong coupling constant, as, 
at HERA together the results of two global analyses. 



a; ~ 0.9. The sensitivity of these cross sections to the gluon was demonstrated 
by evaluating a formed between the measurements and the inclusive jet 
cross section evaluated using various standard PDF sets. The evaluation 
took account of correlations between uncertainties in the variables ij and Et- 
The results indicate that the data prefer a slightly harder gluon distribution 
than is contained within the standard CTEQ4|£3 parameterisation. NLO QCD 
fits incorporating the DO inclusive jet data are being prepared by the MRST 
and the CTEQ groups. The data will provide a stringent constraint on the 
gluon density at high- a;. 

4-3 Structure functions at low-Q"^ 

It is well established that F^^ rises rapidly as x falls below x ^ 0.01 so long 
as is larger than ~ 2 GeV^—apd that NLO QCD gives a good description 
of the data for > 2 GeV^ EJo. Precise measurements of F^'^ for in 
the range 0.045 < < 0.65 GnV^ show that in this kinematic domain 
does not rise rapidly as x fallsEa. The low data can be described using 
non-perturbative models. The x dependence of F^'~^ appears to undergo a 
transition at a of around 1 GeV^ . The nature of the transition between 
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Figure 14. Inclusive double differential single inclusive jet cross sections measured by the 
DO collaboration (points) for jet pseudo-rapidity (r;) intervals plotted as a function of the 
jet transverse energy (Et)- The Standard Model prediction is shown as the solid line. 



the kinematic domain where perturbative QCD is applicable and the domain 
where non-perturbative strong interaction physics is important was discussed 
in a dedicated session at this workshop. The contributions to this session 
are summarised elsewhere in these proceedingsO. A new measurement of 
F^'^ spannirLg.^the range 0.35 < Q"^ < 20 GeV^ was presented by the HI 
collaboratiorO. The low values of were accessed by selecting a sample 
of events in which a photon was radiated in the direction of the incoming 
positron. Such initial state radiation (ISR) causes the effective centre-of-mass 
energy in the e^p collision to be reduced. Hence, a scattered positron of 
particular energy and scattered through a particular angle results from a DIS 
collision is characterised by a smaller value of and a higher value of x than 
would be the case if a scattered positron with the same energy and scattering 
angle were found in the absence of ISR. Figure ^ shows the resulting F2 
plotted as a function of at several values of W, the centre-of-mass energy 
in the virtual photon-proton system. The new HI ISR measurement spans 
the transition region around ^ 1 GeV^. Further measurements in this 
important kinematic region will be required in order that an understanding 
of the underlying physics may be developed. 

The CCFR collaboration has resolved the long-standing discrepancy be- 
tween the structure functions extracted frcpajjacutrino-nucleon CC DIS and 
those extracted from muon- nucleon NC DISO'O. The expression for the cross 
sections for CC vN DIS contains three structure functions, just as Equation 
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Figure 15. Compilation of ZEUS and HI measurements of F2 at low-(3^ in NC DIS (points) 
plotted at several values of W as a function of . The solid lines show the result of the 
HI NLO QCD fit to the HI F2 data for > 3.5 GeV^. The dashed lines show the result 
of the ZEUS Regge fit to the ZEUS F2 data for Q'^ < 0.65 GeV^. 



|lO| does, but the quark composition of these structure functioija^in the Quark 
Parton Model differs 
data were re-analyse^ 



those given in Equations ^ and [l^Zj. The CCFR 



taking charm mass effects into account when eval- 
uating NLO QCD corrections for the longitudinal structure function and the 
difference between the structure function xF^^ for vN and vN CC DIS. In 
addition to an update of the structure function for >|_LpGeV^, the 
CCFR collaboration presented results on F}^ for < 1 GeV^ OlZ3. This is 
the first time that F^ has been determined for such low values of . 

4-4 The longitudinal structure function, F^*^ 

The structure function F^'~^ is multiplied by a factor of in the expression 
for the cross section for NC e~^p DIS. This makes F^'~^ difficult to determine 

at HERA since the NC cross section a'^^ / dxdy oc l/y"^ at fixed x. A direct 
determination of F^^ can be made by reducing the proton beam energy and so 
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raising the value of y for fixed x and 0. The difference between the proton 
beam energies at which HERA has operated to date (820 GeV and 920 GeV) 
is too small to albw an extraction of F^'^ . HI has developed two methods of 
extracting F^'^ EJ. In the first, indirect, method a NLO QCD fit was made 
to the reduced cross section for y < 0.3 where the contribution of F^'~^ may 
be neglected. This fit was then extrapolated to larger values of y and the 
difference between the extrapolation and the data taken as a measure 

0fi^NC_ 

For Q"^ less than 10 GeV^ dn^^/dlny was used instead of the NLO QCD 
fit to extrapolate to high-j/Ea. Secondly, the ISR technique, outlined above, 
has been used to vary the effective centre-of-mass energy so that F^'^ can 
be extracted. Both methods rely upon the measurement of the reduced cross 
section at the highest possible y. At DISOl the HI collaboration presented 
new data on a'^ ^ for 1.5 < < 12 GeV^ from a sample that contains 
scattered positrons with energies as low as 3 GeV. As a result y values as high 
as 0.9 can be accessed. In addition, HI has extracted F^'~^ for values as 
high as 700 GeV^ C3. Figure |l^ shows the values of F^'~^ extracted using these 
methods as well as the F^'^ resulting from the fit to the F2 data for y < 0.3. 
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Figure 16. The longitudinal structure function, Fl, extracted from NC ep DIS data by 
the HI collaboration. 



The structure function R" [F^ w [R^/il + R")] F^) extracted from vN 
and vN data by the CCFR collaboration has been measuredlI3. was pre- 
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sented over the kinematic range 0.004 < a; < 0.5 and 0.5 < < 200 GeY2.. 
The data are well described by NLO QCD evaluated using the MRST PDFsEa. 



. 5 Flavour specific parton density functions 



Much can be inferred about the partonic structure of the proton from the 
inclusive measurements discussed above. A deeper understanding can be de- 
veloped using measurements of the cross sections for processes that are sen- 
sitive to specific quark flavours. The sensitivity of CC data from HERA to 
the d- and u-quark densities was discussed in Section 4.1. An estimate of 



the PDF ratio d/U can be obtained by determining the cross section for the 
Drell-Yan processes pN — > fj,~^ fi~ X where N is either a proton or a deuteron. 
Final measurements of the proton-proton to proton-deuteron Drell-Yaii cross 
section ratio from the E866 experiment at Fermilab were presentedEiJ. The 
Bjorken-x of the two partons that annihilate, xi and X2, can be reconstructed 
from the combined four momentum of the muon pair. When xi » x^ the 
cross section ratio is approximately given by 



d^crPP 



dxidx2 



d{x) 
u{x) 



(15) 



Xi>>X2 



where x — X2- The d/u ratio extracted fmpi Equation 15 by the E866 col- 
laboration is observed to peak at a; ~ O.lSEj. This measurement will give an 
important constraint on the sea quark density at intermediate to high-x. 

Charm production in DIS at HERA has been exploited by both the 
ZEUSe3 and the H1e2I collaborations to extract F^"^, the charm contribution 
to . To date two methods have been used to tag charm in DIS. The 
first method, used by ZEUS and HI, exploits the small mass difference be- 
tween the charmed mesons in the decay D* —> D'^tt to identify D* production. 
The ZEUS collaboration has also used the electrons produced in semi-leptonic 
charm decays. The selection of the semi-leptonic sample was challenging and 
was performed by combining information on the electromagnetic shower pro- 
file in the calorimeter with the energy loss information obtained from the 
central drift chamber. The steps by which was extracted are common to 
both collaborations. First the DIS charm cross section times branching ratio 
was measured in a region of phase space limited by the acceptance of the 
apparatus. Next, the cross section was extrapolated to the full phase space 
using QCD based Monte Carlo programs and the known branching ratios of 
the charmed hadrons taken into account. Finally, .Flf^ was extracted from the 
charm production cross section. The results are plotted in Figure ^ F2'^ 
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is observed to rise rapidly as x falls, a behaviour that mirrors the behaviour 
of the inclusive F^^ . In addition, Fj'^ is observed to exhibit the strong scal- 
ing violations expected in NLO QCD calculations of based on standard 
PDFs. Ff is large, amounting to - 25% of F^'^ for _ 39 GeV^ and 
X ^ 10"'^. The high luminosity that will be provided by the HERA upgrade 
coupled with the increased charm sensitivity offered by the new ZEUS and 
HI silicon micro- vertex detectors promises rapid progress in this field. 




Figure 17. The charm contribution to F^'^ , F^". (a) F^" plotted as a function of x for 
several values of (points). The solid line shows the result of a NLO QCD calculation, 
(b) F2'' plotted as a function of Q-^ for several values of x. The expectations of the HI NLO 
QCD fit are shown as the shaded bands. 



The CCFR and the NuTcV collaborations each presented several new 
results that are sensitive to the flavour content of the proton. NuTeVa pre- 
sented final cross sections for dimuon production in vN and DN DIS. Dimuon 
production in CC i^N DIS is sensitive to the strange quark density. By making 
a NLO QCD fit to the i^N and the PN dimuon data NuTeV has been able to 
demonstrate that the magnitude of the strange sea amounts to approximately 
40% of the up-quark sea and that the s-quark density has approximately the 
same magnitude as the s-quark density. Samples of NC vN and i>N events 
containing a low energy muon were also studied. The muons in these events 
arise from the decay of hadrons produced in the scattering. Single muon pro- 
duction in NC vN and PN DIS is sensitive to the charm quark density. The 
NuTeV results do not require the introduction of intrinsic charm. CCFIO 
has determined the difference between the structure function xF^ in CC DN 
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scattering and the structure function xF!^ in CC vN scattering, A{xF^). In 
the quark parton model /S.{xF!^) = 'ix{s — c). /S.(xF!^) was presented in the 
kinematic range 0.015 < x < 0.08 and for 7 < < 20 GeV^. The data are 
reasonably well described by NLO QCD using standard PDF sets. 

4.. 6 The future of lepton-nucleon deep inelastic scattering 

The upgrade to the HERA accelerator, which was nearing completion at the 
time DISOl was held, will yield a factor five more luminosity and deliver po- 
larised electron and positron beams to the collider experiment ZEUS and HI. 
The ZEUS and HI experiments are also undergoing far reaching upgrades. 
The upgrades themselves and the physics opportunities oilhfi new machine 
and detectors are described elsewhere in these proceedings C3'EJ. Looking fur- 
ther ahead there are two new windows on lepton-nucleon deep inelastic scat- 
tering that were discussed at DISOl. The first, ep collisions between 920 GeV 
protons from the HERA proton ring and 250 GeV electrons or positrons from 
the proposed TESLAJinac, offers exciting possibilities in the study of low- 
X DIS and diffractiorH. The second, vN scattering, would use the intense 
neutrino beam generated by the d£cay in flight of the intense stored muon 
beam of a future Neutrino FactoryEEI. The unprecedented neutrino flux that 
such a facility will provide will allow neutrino targets with a mass of only a 
few kilogrammes to be used. Not only would experiments at such a facility 
be able to determine all the unpolarised neutrino-nucleon structure functions 
with great precision but the flux would be high enough that high precision 
measurements using polarised nuclear targets will be possible. 
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